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1. Introduction 

The MILC collaboration has been carrying out simulations of 2+1 flavor lattice QCD with 
the improved "asqtad" staggered quark action since 1999; for a review of the physics program see 
Ref. [1]. The asqtad ensembles are now complete, and we are progressing to lattice generation with 
the more highly improved HISQ action [2], with promising results so far [3]. Extraction of physics 
with the asqtad ensembles is however continuing. In particular, we are in the process of finishing 
our (asqtad) study of the light pseudoscalar meson sector. Here we give the latest update of this 
project, focusing on SU(3) chiral fits. 

Compared to the last status reports in Refs. [4, 5], we have completed the final four ensembles, 
adding between 35 and 100% more configurations. The ensembles include two of the ones with 
lighter-than-physical strange quark masses, which are crucial for SU(3) chiral fits. This new data 
has allowed us to remove all a priori restrictions on SU(3) NNLO low energy constants (LECs) 
from the fits. It is therefore now possible, for the first time, to report physical values for some of 
these NNLO LECs, although the errors in the values are quite large. In addition, we now consider 
additional alternative versions of the chiral fits, including different treatments of still higher order 
effects. These changes in the analysis give somewhat larger systematic errors on the LECs at NLO, 
as well as on the decay constant and quark condensate in the three-flavor chiral limit. On the other 
hand, errors on fj^ and fx, and on quark masses, are largely insensitive to the changes, and indeed 
have continued to decrease as the ensembles have improved. 

We have also examined in detail the convergence of the SU(3) chiral expansion, as a function 
of quark mass. Results ai^e presented both for the case of three degenerate quark masses and for the 
more physical case of one quark (strange) much heavier than the other two. The expansion appears 
to be significantly better behaved in the former case than in the latter one. 

2. The ensembles and the fitting procedures 

The present analysis uses MILC ensembles at a 0.09 fm, a ^ 0.06 fm and a ^ 0.045 fm. 
Although the MILC collaboration has also generated ensembles at three coarser lattice spacings 
from a 0.18 fm to a 0.12 fm, using that data in the fits is not currently feasible, since it would 
require analytic control of higher order discretization effects. The ensembles used in this study 
are listed in Table 1. The quantity m' is the simulated light {u,d) quark mass. Three ensembles 
have an unphysically light simulated strange quark mass m'^ 0.6/71^ (two a 0.09 ensembles 
with am',. = 0.0186, and one a ?»0.06 ensemble with am'^ = 0.0108), and one ensemble has three 
degenerate quarks with tn' = m'^^ O.lm^ (the a ss0.09 ensemble with am' = am'^ = 0.0031). These 
light-m^. ensembles were created specifically to help control the SU(3) XVY fits. 

As always in MILC staggered simulations, we take the fourth root of the fermion determinant 
in order to eliminate the effect of unwanted degrees of freedom ("tastes") in the lattice generation. 
Recent work (see Ref. [1] and [6] for reviews) makes a strong case that the procedure does indeed 
produce the desired theory in the continuum limit. 

At the fairly small lattice spacings considered here, chiral effects due to taste-violations are 
relatively small but are nevertheless not negligible. We take these effects into account at NLO in our 
chiral fits by using rooted staggered ZPT (rSZPT) [7]. An rSZPT calculation of the two-loop chiral 
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Table 1: List of ensembles used in this study, with mq the tadpole factor and r\/a the scale from the heavy 
quark potential. The ri/a values are "mass-independent," obtained by using an interpolating fit to adjust the 
sea quark masses to their physical values [1]. Errors in these r\ ja values are roughly 3 to 10 in the last digit. 
Configuration numbers with asterisks have been updated since last year 

logarithms does not exist, so at NNLO we use the continuum partially quenched chiral logarithms 
[8], with the the root mean square taste average of the pion mass as argument. In what we call 
"low-mass" fits (using the low-m|^, ensembles, with low valence masses), neglect of taste-violating 
effects is justified at NNLO because such effects are smaller than the mass effects that we keep. 
Thus such NNLO low-mass fits are "systematic" in the sense of ZPT. See Ref. [4] for more details. 

In Ref. [4], the low-mass fits used prior widths to constrain the NNLO LECs to be ^(1) in 
natural units [9] ; LO and NLO LECs were unconstrained. With the current complete data set, we 
can dispense with the constraints and still obtain convergent fits. The data alone now determines 
the NNLO LECs. It is therefore possible to quote physical results for these quantities, although, 
not surprisingly, the eiTors at NNLO are quite large. Eliminating constraints at NNLO also tends to 
increase the errors of the NLO LECs, because there can be more trade-off between the two orders 
in the fit. Clearly the new eiTors are more conservative. Fortunately, eliminating constraints has 
very little effect on quantities extrapolated to physical meson masses, such as the quark masses and 
the decay constants and f^. Such quantities are determined lai^gely by the lattice data (which 
now goes to quite low light quark masses) and are not very sensitive to how the fits divide up 
contributions by orders in ZPT. 

Another feature of the current analysis also tends to increase, relative to Ref. [4] , our estimate 
of errors of certain quantities, in particular, the values of quantities in the 3-flavor chiral limit. 
Define the chiral coupling relevant at two-loops as 1/ (167r^/^^^Q). The value of the decay constant 
Innlo in the coupling is ambiguous a priori because differences in /nnlo (between, say, the 3- 
flavor chiral-limit value, f^, and the physical /;;:) ai^e higher order. In Refs. [9, 1], which did not 
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Figure 1: Convergence of "type-A" NNLO low-mass SU(3) chiral fits for f-^ (left) and mj^/lnix (right) as a 
function of the (degenerate) valence quark mass nix- The strange sea quark mass m'^ is held fixed at w Q.6ms. 
The lines labeled LO, NLO, and NNLO show the ZPT contribution up to and including the indicated order 
The line labeled "+higher orders" shows the effect of adding still higher order analytic terms. The vertical 
line labeled "mass limit" corresponds to the largest value of the valence mass (in the units of the abscissa) 
used in the fits, and is given by Im^ ~ m'^. The fit has been extrapolated to the continuum. 

include NNLO chiral logs, the corresponding issue already arose at NLO. We found there that 
using the "bare" coupling with fyw = /a did not give acceptable fits, and it was necessary to use 
a "physical" NLO coupling (with fn^fNW^fK) to get good fits to our lattice data. Similarly, 
here we cannot obtain good fits of our low-mass data with /nnw = h '■ The best such fits have 
confidence level CL < 0.03, and in addition require unreasonably large lattice-spacing dependence. 
On the other hand, using /nnlo ~ Ik gives good fits, and indeed a value close to fji is chosen by 
the data if /nnlo in the NNLO terms is allowed to be a free parameter. 

Fixing the value of /nnlo is done in two different ways in our fits. Our standard approach 
("type-A fits"), used also in Ref. [4], is to put, for each ensemble, /nnlo = "^/f'' where /|' is 
the decay constant in the chiral limit on that ensemble, and T is a fixed number. The value of T 
is chosen (iteratively) so that /nnlo has the desired value in the continuum limit. For the central 
value fits we ensure f = fj^ m the continuum, and other fits vary /nnlo over a range (roughly 
^■^5fji<fi\iNLO^^-^5fji) that gives acceptable CLs. An alternative approach ("type-B fits") that 
we have tried recently is simply to set, on all ensembles, /nnlo = pfn, where fj^ is taken from 
experiment, and p is a fixed number, again chosen over the range that gives acceptable fits. The 
type-B approach completely decouples the value of /nnlo from the chiral limit quantity that 
appears at lower order and should describe the decay constant in the low-mass regime. 

Chiral fitting is then done in two stages. Systematic low-mass fits are performed first and 
are used to determine low energy constants (LECs) through NNLO: LO parameters /3 and S3 
(often called Fo\/2 and Bq in the literature), NLO Gasser-Leutwyler [10] parameters L,-, and NNLO 
parameters [11, 8] Ki (partially quenched) or Q (unquenched SU(3)). The fits include all partially 
quenched data for pion and kaon (with lighter than physical strange quark mass) decay constants 
and masses (108 points) and have 31 parameters. The full covariance matrix is used, and the CL is 



4 



Results for light pseudoscalar mesons 



C. Bernard 



0.20 



0.15 



— 1 — I — I — I — I — I — I — 

□ a=0.09 fm, alni=0.031 
X a=0.09 fm, amj=0.0186 
K a=0.09 fm, am^^O.OOSl 
♦ a=0.06 fm, am^^O.OlS 
+ a=0.06 fm, amjpO.OlOB 
o a-0.045 fm, am^O.014 

dL = 0.37 
-xV<iof=538/504 (unadj.) 




X 0.0062 
X 0.0031 



K 0.0031 



□ 0.0124 

□ 0.0093 

□ 0.0062 
0.00465 

□ 0.0031 

□ 0.00155 



+ 0.0072 

full, conl., + 0.0036 + 0.0054 

* extrap + 0.0036 

iiiexpt. (r,|=0.3133fm) 0.0026 0.0025 
_i I I I I I i_ I I I t 0.0018 



0.00 



0.05 
2m„r 1 



0.10 

X (Zn,/zr) 



0.15 



7.5 



^ 7.0 

X 

S 
w 

6.5 



6.0 — 



0.00 



— I 1 1 1 1 1 — 

□ 0=0.09 fm. am^O.031 
« a=0.09 fm. amj^O.OOSl 

_ X a=0.09 fm. am|i=0.0186 
+ a=0.06 fm, amJ^O.Ola 
+ a=0.06 fm. am^O.0108 
o a=0.O45 fm, am'^0.014 

dL = 0.37 

//dof=538/504 (unadj.) 



X 0.0068 
X 0.0031 



□ 0.0124 

□ 0.0093~ 

□ 0.006E - 
0.00485 

□ 0.0031 " 

□ 0.00155- 

♦ 0.0072 
+ 0.0054 

♦ 0.0038- 
0.0025 

♦ 0.0018 

o 0.0028 




0.05 0.10 
2m,r, X (Z^/zr) 



Figure 2: High-mass, type- A, SU(3) chiral fits, using ri = 0.3133. We show partially quenched data points 
with degenerate valence masses for the decay constant (left) and for m^/ [Inix) (right). The red lines show 
the results after extrapolation to the continuum, putting m[ — nis, and setting light sea and valence masses 
equal. The decay constant is also extrapolated to the physical mass point and compared with experiment 
(blue burst). 



good (0.50). The convergence of XPT in these fits (type-A) is shown in Fig. 1. Although the NLO 
terms for m^/lnix seem to be anomalously small, the overall convergence is satisfactory. 

In the second stage of fitting, the "high-mass" fits, all ensembles listed in Table 1 are included 
with the valence masses restricted to nix + my < 1. 2ms. In the type-A case, the LO, NLO, NNLO 
LECs are fixed (for central values) or fixed within statistical errors at the values from the low-mass 
fits. N^LO and N'^LO analytic terms are included, but not the corresponding logs. These terms are 
needed to obtain good CLs, and they allow us to interpolate around the (physical) strange quark 
mass. The fact that they are required indicates that SU(3) XPT is not converging well at these 
mass values, unlike the situation in the low-mass case (see below for more evidence). However, 
since the LO, NLO, and NNLO LECs dominate the chiral extrapolation to the physical point, the 
results for decay constants and masses are rather insensitive to the form of these N^LO and N^LO 
interpolating terms, as long as the fits are good. 

In the type-B case, high-mass fits with the LO, NLO, and NNLO LECs fixed at their low-mass 
values do not give good CLs. We can obtain acceptable fits by allowing these LECs to vary, with 
prior widths set by the low-mass statistical errors. But in that case, several of the LECs move by 
more than 2a. For these reasons the type-B fits are currently disfavored, although we do include 
them as alternative chiral fits in estimating systematic errors. 

The (central value, type-A) high-mass fits are illustrated for pions (degenerate valence points) 
in Fig. 2, where the scale has been set by ri = 0.3133(23) fm [12]. The extrapolations to the 
continuum and to the physical masses are also shown. The agreement of our result, = 129.2 it 
0.4 lb 1 .4 MeV, with the experimental value (130.4(2) MeV [13]) is good. From now on we set the 
scale using fji, which gives ri = 0.3 106(8) (14) (4) fm, where the last error is experimental. 

Figure 3 shows the convergence of XVY for the pion decay constant. The left plot is for the 
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Figure 3: Convergence of ZPT for the decay constant, with the same fit as in Fig. 2. Left, the degenerate 
case (m^. = m' = nix) as a function of m^. Right, 2-flavor chiral hmit nix = m' = &s a function of m'^. These 
fits have been extrapolated to the continuum. 



degenerate 3-flavor case. Note that the convergence is still good significantly beyond the point 
where 2mx = rus, where the meson mass is already close to the physical value of tuk- The right 
plot gives the decay constant as a function of the strange sea quark mass m'^, with the valence and 
light sea masses extrapolated to zero (the 2-flavor chiral limit). The vertical line indicates = m^, 
when the kaon is near physical. Here the convergence is reasonable only up to about m^. = 0.6m.s; 
by m'^ = 0.8m.s the NNLO term has the wrong sign, and moves the NLO result away from the 
"higher orders" line (which must be close to what the lattice data demands, since the fit is good). 
Nevertheless, the NLO contribution remains reasonable; the problem is mainly at NNLO. This kind 
of behavior is not unexpected for an asymptotic expansion. 

For the m\/{2mx) case (not shown), the convergence of ZPT is similar to that of the decay 
constant. Again, convergence in the degenerate 3-flavor case is good up to meson masses about 
10% or 15% higher than the kaon. In the 2-flavor chiral limit, the behavior as a function of m[ 
appears to be comparable to that of the decay constant, breaking down before m[ = mg. 



3. Preliminary results 

The following is a sample of our current results: 

= 156.1 ±0.4 +ojMeV, /^/^ = 1.197 (2) (+3) , 

K, = 0.2247 (+'^) , /3 = 113.6±3.6±7.7MeV , 

/2 = 123.0 ± 0.5 ± 0.7 MeV, 2U- Ls = -0.5\{n){t%) , 

2L6-L4 = 0.09(24) (+32) , Lj = 1.79(16)(+28) , 

L4 = 0.l9{22){tll) , 2C21 -Ci9 = 1.5(6)(1^) , 

Kig = 3.5 lb 1.2 '^^j , A'39 —Kn = 3.4 it 1.5 +};g , 
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h = 3.18(50)(89) , /4 =4.29(21)(82) , 

where errors are statistical and systematic. /2 is the decay constant in the 2-flavor chiral limit. The 
NLO LECs Li are in units of 10"^ and the NNLO LECs Ki and C, are in units of 10"^; both are 
at chiral scale m,j . Other Ki and C, are also of this order of magnitude, but most have statistical or 
systematic errors that are more than 100%. The scale invariant SU(2) LECs £3,4 are obtained from 
the SU(3) LECs using the two-loop conversion formulae [14]. There is good agreement between 
the SU(3) chiral fit results described here and the results of the SU(2) chiral fits [5, 15] for all 
quantities that can be directly compared. 

As discussed above, errors for /3 and the NLO LECs are larger than previous ones [4] , while 
those for fx, Ik/ fn, fi, as well as for quark masses (not given here), are similar to or somewhat 
smaller than before. In general, we use type-A fits with /nnlo — fn for central values, and alterna- 
tive chiral fits (including type-B) for systematic error estimates. The exception is f-i, where we have 
averaged type-A and type-B results and symmetrized the eiTors. This is because the distribution 
for /3 from type-A and type-B fits is bimodal, with the former giving values toward the high end 
of the range, and the latter giving values toward the low end. The large uncertainty in the 3-flavor 
chiral limit, compared to the physical point, comes ultimately from the fact that we only have one 
usable ensemble with three degenerate light flavors, while there are many with m[ 171^. 

We thank J. Bijnens for his program to compute the partially quenched NNLO chiral logs. 
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